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Abstract 

We investigate the singularity formation of a 3D model that was recently proposed 
by Hou and Lei in [8] for axisymmctric 3D incompressible Navier-Stokes equations with 
swirl. The main difference between the 3D model of Hou and Lei and the reformulated 
3D Navier-Stokes equations is that the convection term is neglected in the 3D model. 
This model shares many properties of the 3D incompressible Navier-Stokes equations. 
One of the main results of this paper is that we prove rigorously the finite time sin- 
gularity formation of the 3D model for a class of initial boundary value problems with 
smooth initial data of finite energy. We also prove the global regularity for a class of 
smooth initial data. 

Key words: Finite time singularities, nonlinear nonlocal system, stabilizing effect of 
convection, incompressible Navier-Stokes equations. 

1 Introduction 

The question of whether a solution of the 3D incompressible Navier-Stokes equations can 
develop a finite time singularity from smooth initial data with finite energy is one of the 
most outstanding mathematical open problems [H [T3l I16j . Most regularity analysis for 
the 3D Navier-Stokes equations relies on energy estimates. Due to the incompressibility 
condition, the convection term does not contribute to the energy norm of the velocity field 
or any (1 < p < oo) norm of the vorticity field. As a result, the main effort has been to 
use the diffusion term to control the nonlinear vortex stretching term without making use 
of the convection term explicitly. 
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In a recent paper by Hon and Lei [8], the authors investigated the stabiUzing effect of 
convection by constructing a new 3D model for axisymmetric 3D incompressible Navier- 
Stokes equations with swirl. Specifically, their 3D model is given below. 



dtui = u{dl + -dr + dl)ui + 2dztpiui, 
r 

dtui = v{dl + -dr + dl)uji + d,{{uif) 

r 

-{dl + -dr + dl)i^i=uji. 



(1) 



(2) 



(3) 



Note that ([I])-© is already a closed system. The only difference between this 3D model 
and the reformulated Navier-Stokes equations is that the convection term is neglected in the 
model. If one adds the convection term back to the left hand side of ([T]) and one would 
recover the full Navier-Stokes equations. This model preserves almost all the properties of 
the full 3D Navier-Stokes equations, including the energy identity for strong solutions of 
the 3D model and the divergence free property of the reconstructed 3D velocity field given 
by = rui, = — 52(r-0i), = }:dr{r'^4'i)- Moreover, they proved the corresponding 
non-blowup criterion of Beale-Kato-Majda type [1] as well as a non-blowup criterion of 
Prodi-Serrin type [HI [15] for the model. In a subsequent paper, they proved a new partial 
regularity result for the model [9] which is an analogue of the Caffarelli-Kohn-Nirenberg 
theory [2] for the full Navier-Stokes equations. 

Despite the striking similarity at the theoretical level between the 3D model and the 
Navier-Stokes equations, the former has a completely different behavior from the full Navier- 
Stokes equations. In [8], the authors presented numerical evidence which supports that 
the 3D model may develop a potential finite time singularity. They further studied the 
mechanism that leads to these singular events in the 3D model and how the convection 
term in the full Navier-Stokes equations destroys such a mechanism (see also [5]). 

One of the main results of this paper is that we prove rigorously the finite time singularity 
formation of this 3D model for a large class of initial boundary value problems with smooth 
initial data of finite energy. In our analysis, we consider the initial boundary value problem 
of the generalized 3D model which has the following form [8] (we drop the subscript 1 and 
substitute ^ into 



where A is a n-dimensional Laplace operator with (x, z) = (xi, X2, 2:). Our results 
in this paper apply to any dimension greater than or equal to two {n>2). To simplify our 
presentation, we only present our analysis for n = 3. We consider the generalized 3D model 
in both a bounded domain and in a semi-infinite domain with mixed Dirichlet and Neumann 
boundary conditions. One of the main results of this paper is the following theorem. 




(4) 
(5) 
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Theorem 1.1 Let = (0, a) x (0, a), 17 = x (0, b) and T = {(x, z) \ x £ Qy_, z = 0}. 
Assume that uq > for (x, z) G 17, uo\gQ = 0, and uo,xpoz G H'^{^). Moreover, we assume 
that ip satisfies the following mixed Dirichlet Neumann boundary condition: 

V'lan\r = 0, (^, + /3V')|r = 0, (6) 
u < a < \/'ZT:/a ana 'i 

condition: 



^^ih /3 > ^ (f^f^ j ■ Define y, .) = J sin (^) sin (^) «;/.ere 

a satisfies < a < \/27r/a and 2 (^)^ J^^^'^+e-"'') ~ f^' V uq and ipo satisfy the following 



/ (loguo)^(ix(iz > 0, / ipoz<pdxdz > 0, (7) 
in Jn 

i/ien the solution of the 3D model (^-(0^ will develop a finite time singularity in the 
norm. 

The analysis of the finite time singularity for the 3D model is rather subtle. The main 
technical difficulty is that this is a multi-dimensional nonlinear nonlocal system. Currently, 
there is no systematic method of analysis to study singularity formation of a nonlinear 
nonlocal system. The key issue is under what condition the solution u has a strong alignment 
with the solution ip^ dynamically. If u and tp^ have a strong alignment for long enough 
time, then the right hand side of the u-equation would develop a quadratic nonlinearity 
dynamically, which would lead to a finite time blowup. Note that is coupled to u in 
a nonlinear and nonlocal fashion. It is not clear whether u and ipz will develop such a 
nonlinear alignment dynamically. As a matter of fact, not all initial boundary conditions 
of the 3D model would lead to finite time blowup. One of the interesting results we obtain 
in this paper is that we prove the global regularity of the 3D inviscid model for a class of 
smooth initial boundary data. Moreover, we obtain a relatively sharp global bound on the 
H^-noim (s > 4) of the solution in terms of its initial condition. 

One of the main contributions of this paper is that we introduce an effective method of 
analysis to study singularity formation of this nonlinear nonlocal multi-dimensional system. 
There are several important steps in our analysis. The first one is that we renormalize the 
u-equation so that the right hand side of the renormalized u-equation becomes linear. This 
is accomplished by dividing both sides of dH) by u and introducing log(u) as a new variable. 
This is possible since tig > in 17 implies that u > in 17 as long as the solution remains 
smooth. The renormalized system now has the form: 

{log{u))^ = 2^z, (x,z)Gl7, (8) 
-AVt = (n^)^. (9) 

This idea is similar in spirit to the renormalized Boltzmann equation introduced by DiPerna 
and Lions in their study of the global renormalized weak solution of the Boltzmann equa- 
tions |3]. The second step is to work with the weak formulation of the renormalized model 
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([H])-® by introducing an appropriately chosen weight function 4> as our test function. How 
to choose this weight function (p is crucial in obtaining the nonlinear estimate that is re- 
quired to prove finite time blowup of the nonlocal system. Guided by our analysis, we look 
for a smooth and positive eigen-function in O that satisfies the following two conditions 
simultaneously: 

- = Ai(?!>, d'^(f) = X2(l), for some Ai, A2 > 0, {x,z)eQ. (10) 

The function (p defined in Theorem 1.1 satisfies both of these conditions. We remark that 
such eigen-function exists only for space dimension greater than or equal to two. In the 
third step, we multiply to ([8]) and (j)z to ([9]), integrate over and perform integration by 
parts. We obtain by using (fTOl) that 

d 



dt 



/ (logu)0dxdz = 2 ipz(j)dxdz, (11) 
Jn Jn 



d 



Ai^ / V'^(/)dxdz = A2 / u^^dxdz. (12) 
dt Jn Jn 

All the boundary terms resulting from integration by parts vanish by using the boundary 

condition of ip, the fact that u\z=o = u\z=b = 0, the property of our eigen-function cj), and 

the specific choice of a defined in Theorem 1.1. Substituting into (fTTj) gives the crucial 

estimate for our blowup analysis: 

d^ f 2A2 

Further, we note that 



/ (log n)0dxd2; = f u^(j)dxdz. (13) 
Jn Ai Jn 



log(ti)(/>dxd2; < / (log(n))'^(/)dxdz < / u^dxdz 
Jn Jn 



\ 1/2 / ^ \ 1/2 2a / f \ ^^"^ 

< ( I (/.dxdzj i^J (/.M^dxdzJ = -—j= i^J ^u^dxdzj . (14) 



Integrating (jl3p twice in time and using (jl4p . we establish a sharp nonlinear dynamic 
estimate for (J^ (pu'^dxdz)^^'^ , which enables us to prove finite time blowup of the 3D model. 

This method of analysis is quite robust and captures very well the nonlinear interaction 
of the multi-dimensional nonlocal system. As a result, it provides a very effective method 
to analyze the finite time blowup of the 3D model and gives a relatively sharp blowup 
condition on the initial and boundary values for the 3D model. Using this approach, we 
prove finite time blowup of the 3D model for a large class of initial boundary value problems 
in both bounded and semi-infinite domains. We have extended this method of analysis to 
study a generalized 3D model by changing the sign in front of the Laplace operator in ([5]). 

We also study singularity formation of the 3D model with /3 = in ([6]). This case is 
interesting because the strong solution of the corresponding 3D model satisfies an energy 
identity. In this case, we can establish a finite time blowup under an additional condition : 

j j {ip - ipo)\r sin (J^^ sin (J^^ dx < Co J Voz^dxdz, 
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as long as the solution remains regular, where cq > depends only on the size of the domain. 

We remark that although the 3D model using the mixed Dirichlet Neumann boundary 
condition with (3^0 does not conserve energy exactly, we prove that the energy remains 
bounded as long as the solution is smooth and (5 < cq for some cq > 0. We also establish 
the local well-posedness of the initial boundary problem with the mixed Dirichlet Neumann 
boundary condition. Our numerical study shows that the energy is still bounded up to 
the blowup time even if /? > cq. Our numerical study also suggests that the nature of the 
singularity in the case of /? > co is qualitatively similar to that in the case of /? < cq. 

Another interesting result is that we prove finite time blowup of the 3D model with 
partial viscosity. Under similar assumptions on mq) V'o and ojq as in the inviscid case and by 
assuming that oo satisfies a boundary condition similar to ^Z', we prove that the 3D model 
with partial viscosity 

Ut = 2uipz, (15) 
uJt = {u^)^ + iyAoj, (16) 
-Ai/j = uj, (17) 

develops a finite time singularity. We conjecture that the fully viscous 3D model would 
also develop a finite time singularity. For the fully viscous model, we need to propose a 
boundary condition for u. Currently we have not been able to prove finite time blowup of 
the fully viscous model using our method of analysis unless we impose a Dirichlet-Neumann 
boundary condition on the renormalized variable log(ii), which is not very natural. 

We remark that Hou, Li, Shi, Wang and Yu [lOj have recently proved the finite time 
singularity of a one-dimensional nonlinear nonlocal system: 

Ut = 2uv, vt = H{u^), (18) 

where H is the Hilbert transform. This is a simplified system of the original 3D model 
along the symmetry axis. Here v plays the same role as ipz- The singularity of this nonocal 
system is remarkably similar to that of the 3D model. The stabilizing effect of convection 
has also been studied by Hou and Li in a recent paper [7] via a new ID model which can 
be used to construct a family of exact solutions of the 3D Euler or Navier-Stokes equations. 
They showed that convection plays an essential role in canceling the destabilizing vortex 
stretching in the new ID model and obtained an a priori pointwise estimate for a Lyapunov 
function which characterizes the derivatives of the vorticity in their model. 

The rest of the paper is organized as follows. In Section 2, we study some properties of 
the 3D model. In section 3, we prove the finite time blowup of the inviscid 3D model with 
mixed Dirichlet Neumann boundary conditions. In Section 4, we analyze the finite time 
blowup of the 3D model with partial viscosity. Section 5 is devoted to analyzing the finite 
time blowup of the 3D model with boundary conditions that conserve energy. We prove the 
global regularity of the 3D model for a class of smooth initial boundary value problems in 
Section 6. A technical lemma is proved in Appendix A. 
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2 Properties of the 3D model 

2.1 Local well-posedness in H'^ 

In this section, we will establish the local well-posedness of the initial boundary problem of 
the 3D model with the mixed Dirichlet Neumann boundary condition. We will present our 
analysis for the semi-infinite domain using the Sobolev space H'^. The same result is also 
true in a bounded domain. 

Consider the 3D model with the following mixed initial boundary condition: 

I ut = 2u^, (x,z)eO = J7,x(0,oo) (19) 

V'|t=o = V'o(x, 2;), u\t=o = uo{x,z) >0, 

where x = (xi, X2), J^x = (0, a) x (0, a), F = {(x, z) | x G ilx, z = 0} and A = Ax + -q^ = 

92 92 92 

The local well-posedness analysis depends on an important property of the elliptic op- 
erator with the mixed Dirichlet Neumann boundary condition. 

Lemma 2.1 There exists a unique solution v G H^{^}) to the boundary value problem: 

-Av = f, (x,z)gJ7, (20) 
\an\r = 0, {v, + Pv)\r = 0, (21) 



V 



if P £ Soo ^ {P \ P ^ ^ for all k £ 1?}, f G H'-'^{n) with s > 2 and f\dn\r = 0. 
Moreover we have 

\\v\\H^in)<Cs\\f\\Hs-2^n), (22) 
where Cg is a constant depending on s, \k\ = \/k'f + 1^2- 
We defer the proof of Lemma 12.11 to Appendix A. 

Remark 2.1 We remark that we can prove the same result as in Lemma \2.1\ for a bounded 
domain VL = ilx x (0, b) with the same boundary condition by assuming that P £ Sb where 

Tr\k\ Tr\k\ / 1 -I- p-2|*:kV«\ 

S, = {P\P^Mandp^^i j^^^^^ for allk £ I?}. (23) 



We also need the following well-known Sobolev inequality |12j . 
Lemma 2.2 Let u,v £ H'^{i}) with s > 3/2. We have 

WuvllH^in) < c\MH'={n)\\v\\H^n)- (24) 
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Now we can state the local well-posedness result for the 3D model with the mixed 
Dirichlet Neumann boundary condition. 

Theorem 2.1 Assume that no,'i/'Oz £ s > 3/2, and uolan = 0. Moreover, we 

assume that (3 G 6*00 (or Sb) as defined in Lemma \2.1\ and satisfies the following mixed 
Dirichlet Neumann boundary condition: 

i^\dn\r = 0, (V.+/3V')|r = 0. (25) 

Then there exists a finite time T = T (||uo||//s(n); ||V'0z||_f/s(f2)) > such that the system 
with boundary condition I125\} has a unique solution, Ujtpz S C^([0, T), 77^(17)). 

Proof Let w = , then we obtain an equivalent system for w and ijj as follows: 

wt = 4wij,, (26) 
-Ailjt = w,. (27) 

Using Lemma 12.11 and Lemma 12.21 we get 

\\wt\\H'> < Cs||'u;||//.||^/;^||h'', (28) 

\\lpt\\H^+^ < CsWWzWh^-'' < Cs\\w\\h''- (29) 

On the other hand, for any a = {ai, 02, 03) with aj > (j = 1, 2, 3), we have 

= 2(9"z;t,5%) <2||a%||i2||a'^T;||i2. (30) 



Summing over |a| < s gives 

^ '\v\\hs < 2\\vt\\H'>\\v\\H'', (31) 



dt' 



which implies 



dt 



Using ([52|) . we obtain 



Substituting §3^ into (128])- ([29]), we get 

< CsII'wIIh-IIV'zIIh'', (34) 

^IIV'zIIh^ < Cs\\w\\h'^- (35) 

From the above a priori estimate, we can prove the local well-posedness of the 3D model 
with boundary condition (I25p . □ 

Remark 2.2 We remark that the above analysis can be easily modified to prove the local 
well-posedness of the initial boundary value problem of the partial viscous system I^15\)-{n) 
for initial condition uq, V'Oz G H^{Q) and uq S H^~^{Q) with s > 3/2 (see also fB^). 
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2.2 Bounded energy for the 3D model with mixed boundary conditions 

Proposition 2.1 Let Q. = ilx x (0,^) and T = {(x, z) | x G ilx; z = 0}. Assume 
uq\z=q = UQ\z=b = 0, uo,ipoz G H^{i^) andil^ satisfies the following mixed Dirichlet Neumann 
boundary condition: 

^\dn\r = 0, (V.+/?V')|r = 0. (36) 
Then the strong solution of the 3D model satisfies the following identity 



d 
It 

Moreover, we have 



[ {u^ + 2\Vi;\^) d^dz - 2(3 [ [ V'^|^=orfx ) = 0. (37) 
Jn Jo Jo J 



[ {u^ + 2{1- pb)\Vi;\^)d^dz < [ {ul + 2\V'4>o\'^) di^dz - 2(3 f [ VoU=odx. (38) 
Jn Jn Jo Jo 

Remark 2.3 One immediate consequence of the above proposition is that if [3 < 1/6, both 
f^u'^d'Kdz and | V^'pc^xdz are bounded. 

Proof Multiplying by n and integrating over yields 

^ / n^dxdz = 4 / u^V'^dxdz. (39) 



n 



dt 

Next, we multiply ([5]) by and integrate over 17 to obtain 

- / AV'tV'dxdz= / (u^)^Vdxdz. (40) 
Jn Jn ^ 

Integrating by parts and using boundary condition (j25p . we have 

/ |VV'pdxdz + 2 / 'i(jzt'4^\z=odx = -2 / uV^dxdz. 
Jn Jn^ Jn 

Multiplying (|^T|) by 2 and adding the resulting equation to gives 



(41) 



^ [ (u2 + 2|VVp)dxdz 
dt Jn 



dt 

which gives (I37p . On the other hand, we have the following estimate 

Jn^ \Jo 
< b 



4/3 / ^j^|^=odx 
Jn^ 

2/3^ / V''|.=odx, (42) 



V'^|^=odx = / I / t(j^dz ] dx 



> [ [ V^dzdx < b [ iViZ-pdxdz. (43) 
Jn^ Jo Jn 



This implies that 

[ (u2 + 2|VV'P)dxdz-2/? /" ^2|^=odx> [ (u^ + 2(1 - V^/;^) dxdz. (44) 
Jn Jn^ Jn 

Combining (j37p with ()44p . we obtain 

/" ('u2 + 2(l-/36)|V^|2)dxd2< / (ug + 2|V^oP)dxdz-2/? /" ^oU=odx. (45) 
This completes the proof of Proposition 12.11 □ 

3 Blow-up of the 3D inviscid model 



In this section, we will prove that the 3D model (jl9p develops a finite time singularity for 
a class of smooth initial data with finite energy. The finite time blowup is proved in a 
semi-infinite and a bounded domain with mixed Dirichlet-Neumann boundary conditions. 

3.1 Blow-up in a semi-infinite domain 

First, we consider the initial boundary value problem (|19p-()20p in a semi- infinite domain 
with 17 = X (0, oo). The main result is stated in the theorem below: 

Theorem 3.1 Assume that uq\qq = 0, uq\q > and UQ,ipQz £ H'^(Q,). Further we assume 



V2 



TT 



that P > and (5 G 5*00 as defined in Lemma [2A[ Choose a = and define 

a 

(j){x, z) = e~°'^(j)i{x.), (/)i(x) = sin sin , {x,z)£^}, (46) 

a a 

f f VTQ^/^ f°° dx 

A= (logMo)'Adxdz, B = 2 ipo^cpdyidz, D = — , loc = , „ ■ 

Jn Jn a(2(f) -a2) Jo Vx^ + l 

If A > and B > 0, then the 3D model il9\) -^2U \) will develop a finite time singularity in 
the H^-norm no later than 

(2DB ./3tTT\~^'^ 

-'oo • 



I- 



2a 



Proof By Theorem 12. H we know that there exists a finite time T > such that the 
system (fT9]) has a unique smooth solution, u^ipz £ C'^([0, T), //^(ri)). Let T^, be the largest 
time such that the system (|19p with initial condition uq,iPq has a smooth solution m, V'z £ 
C^{[0,Tf,); H'^{i^)). We claim that Tf, < oo. We prove this by contradiction. 

Suppose that Th = oo, this means that for the given initial data uo,tpo, the system (fT9|) 
has a globally smooth solution u, ipz £ C^([0, oo); H'^{^)). Multiplying (pz to the both sides 
of ([5j) and integrating over we get 

A^Pt(pzdyidz = / (u^)6zdyidz. (47) 



n 
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Note that u\Qn = as long as the solution remains smooth. By integrating by parts and 
using the boundary condition on ip and the property of (j) to eliminate the boundary terms, 
we have 

- / V'^t A^dxdz - / 'il^zt(f>z\z=odx- / ViAx0U=odx= / n^</)^^dxdz. (48) 
Substituting (p into the above equation, we obtain 

.2 



= a" u^d^dz + /^^ (q/3 - ^) ^t|,=o0i(x)dx. (49) 



(/)i(x)dx 

2=0 



By the definition of a, we have 

27r2 ^ , 27r2 ^/27r 
= 0, and a = — 

pa 



a/3 ^ = 0, and a = -—7 < , (50) 



since (3 > Thus the boundary term on the right hand side of ([49l) vanishes. We get 



— / ^pz(pdxaz = rr 



— / Vz^dxd^ = — / ^^<?!'dxd2. (51) 



d_ 

dt 



Next, we multiply (/> to ([8|) and integrate over fi. We obtain 

■ / (logu)</>dxdz = 2 / V^^dxdz. (52) 
Combining (|5ip with (j52p . we have 

— 77 / (log M)(?!)dxd2; = / u (T^dxdz. (53) 

dt^Jn^ 2(f)'-a2 7Q ^ ^ 

Integrating the above equation twice in time, we get 

/ (logu)(/)dxdz = f f ( f u'^(l)d:x.dz^ drds + A + Bt 

Jn 2 {^f - q2 Jo Jo \Jn J 

> [ [ ( f u^</'dxdz^ drds + Bt. (54) 

- 2(f)' -02 70 Jo \Jn J ^ ' 

Note that u > for (x, z) E and t < Th. It is easy to show that 

/ (logu)(pd'Kdz < / (log u)^</>dxdz < / n(^dxdz 
Jn Jn 

1/2 / . ^ 1/2 

< ( / (*dxdz 1 / n^Mxdz 



2a 



1/2 

u'^4>dxdz , (55) 



10 



where {logu)~^ = max(logti, 0). Combining ()54p with ()55p gives us the crucial nonhnear 
dynamic estimate: 



Define 



^(*) = ,^ !^^J [ [ ( L ^^Vdxdz") drds + ^^^Bt. (57) 



a(2(f) -a2) Jo Jo VJn 



/ TT 

Then we have -F(O) = and = > 0. By diff'erentiating (1570 twice in time and 

2a 



substituting the resulting equation into (j56p . we obtain 



dt2 



/ n^dxdz > DF"^, (58) 
in 



where = ^ . Note that Ft = D {j^ u^^dxdz) ds + ^B > 0. Multiplying 

a(2(f) -q2) 

Ft to and integrating in time, we get 



2 

where C = (Ft(0)f = ^B^. Define 



Then, integrating (|59p in time gives 



(60) 



Note that both / and F are strictly increasing functions, and I(x) is uniformly bounded for 
all > while the right hand side increases linearly in time. It follows from (j60p that F{t) 
must blow up no later than 



J 



T — T 



This contradicts with the assumption that the 3D model has a globally smooth solution. 
This contradiction implies that the solution of the system (I19p must develop a finite time 
singularity no later than T* . □ 
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Remark 3.1 As we can see in the proof of Theorem \3.1l the same conclusion still holds if 
we replace the boundary condition 

(V'. + /3^)U=o = 0, 

by the following integral constraint 

/ [Wz + P'>P)\z=oSin sm dx = 0. 

3.2 Blow-up in a bounded domain 

In this subsection, we will prove finite time blow-up of the 3D model in a bounded domain. 
First, we formulate the initial boundary problem of the 3D model as follows: 

{-a2 = (n't ' (^'^)^" = ^-^(0'^)' (6^) 
i^\dn\r = 0, {^Pz + Pi^)\r = 0, (62) 
V'|t=o = ^o(x, z), u\t=o = tio(x, z) > 0, 

where x = (xi,X2), f^x = (0, a) x (0,a), T = {(x, z) G 17 | x G $7x; z = 0}. We can get a 
similar blow-up result which is summarized below: 

Theorem 3.2 Assume thatuQ\QQ = 0, uq\q > 0, and uo,tpoz £ if^(ri). Further, we assume 
that P G Sb as defined in Lemma \2.1\ and satisfies (3 > — p P . Define 

a \ g\/27rf)/a g— V27rfe/a / 

(p(x, zj = sm sm , (x, (b3) 

2 a a 

where a satisfies < a < \P1t\ ja and 2 (^)' Q,(e°'M^e-°'') ~ ^' 

5/2 .oo 



f f ■na^l^ f 

A= (logno)(Adxdz, B = 2 Voz^dxdz, D= —p, , loo = 

Jn Jn a(2(f) -a2) Jo Vx^ + l 

If A > and B > 0, then the solution of ^61l)-^6^ will blow up no later than 



V 3 2a 



1/3 

'oo ■ 



Proof We follow the same strategy as in the proof of Theorem 13. 1[ By Theorem 12. H 
we know that there exists a finite time T > such that the system (j6ip has a unique 
smooth solution, u,^pz £ C^{[0,T), H^{Q)) for < t < T. Let Tf, be the largest time 
time such that the system (I19p with initial condition UQ,ipQ has a smooth solution u, V'z S 
C-^{[0,Tb); H'^{^)). We claim that Tb < oo. We prove this by contradiction. 
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Suppose that = oo, this means that for the given initial data uo,tpQ, the system ([6T]) 
has a globahy smooth solution u, ipz S C^([0, oo); H^{Q,)). Multiplying (pz to the both sides 
of the V'-equation and integrating over J7, we get 

- / A^JtcPzdyidz = I {u^)^(Pzdxdz. (64) 
Jn Jn 

Note that u\Qn = as long as the solution remains smooth. By integrating by parts and 
using the boundary condition on ^ and the property of (p to eliminate the boundary terms, 
we have 

- / V'^tA^dxdz - / 'il^zt(f>z\z=odx- / VtAx0U=odx= / n^</)^^dxdz. (65) 
Substituting </> to ()65p and using the boundary condition for ip, we obtain 



(Ai(x)dx 

z=0 



/ n^^dxdz + / ( ^ fe"'' - e""'') /? - ^ fe"' + e""^ ) ) V'tU=o'/'i(x)dx 



2 V ; ^ a2 



/ nVdxdz + ^ fe-" - e^"") / ('^ " ^ (3^ - ""l^'^rl ) ^tU=o0i(x)dx, (66) 



2V ' Jn^\ ^a' a(e"''-e-°*) 



,2 eaf- _)_ e~°^'' 

where = sin sin Let h(a) = 2 ( — j — j—^ -ab'j ' ^^^^^^ computations 



show that ■^h{a) < for all a > 0. Thus we have 

^ ^— \=h \ < h{a) < /i 0+ = cx), < a < . (67) 

a \QV2nb/a _ ^-V2nb/a J \ a J ^ ' ^ a ^ ' 

Since P > h (^^) by assumption, we can choose a unique a with < a < such that 

/?. (68) 



27r2 e"'' + e""'' 



a2 Q!(e"^ - e-"^y 

With this choice of a, the boundary term in (j66p vanishes. Therefore we obtain 

d /" (x^ f 

— / V-'z^dxdz = 2 / ^i^^dxdz. (69) 



dt Jn"^" 2(^)^-a2 Jn 

Next, we multiply to ([5]) and integrate over 0. We get 

■ / (logu)</>dxdz = 2 / ipzcpdxdz. (70) 
Jn 



d 
di 
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Combining (j70|) with (|69)) . we get 

f 2c? f 

^y^(logn)0dxd. = -^^^y^nVdxd.. (71) 

Now we can follow the exactly same procedure as in the proof of Theorem 13.11 to prove that 
the 3D model must develop a finite time blow-up. □ 

Remark 3.2 We remark that the same conclusion is still true if we replace the Dirichlet 
boundary condition '4^\z=b = by the Neumann boundary condition ■4'z\z=b = 0. The only 
difference is that the weight function (p is now changed to 

(pix,zj = sm sm , (x, zlGSZ, (72) 

2 a a 

where < a < \pl-K ja, and (3 satisfies a variant of I123\) in Lemma \2.1\ and 

' * < /5 < 26 - . (73) 



a \ ^V^ixb/a _|_ g-V27rfe/a / \ a 

We omit the proof here. 

3.3 Blow-up of a generalized 3D model 

In this section, we study singularity formation of a generalized 3D model by changing the 
sign of the Laplace operator in the •^-equation ([S]). Specifically, we consider the following 
generalized 3D model: 



= iu 



2s , (x,z) e = X (0,a) = (0,a) X (0,a) X (0,a). (74) 



The boundary and initial conditions are below 

^l90\r = 0, -iA^lr = 0, r = {(x, z) | x G fi^, z = Q, oi z = a} (75) 
ip\t=o = ipQ{yi,z), u|t=o = no(x, z) > 0. 

In this subsection, we will generalize the singularity analysis presented in the previous 
subsection to prove that the solution of the generalized 3D model will develop a finite time 
singularity. The main result is summarized in the following theorem. 

Theorem 3.3 Assume that uq\q^ = 0, uq\q, > 0, and uq^iJjqz S H'^{^). Further, we define 

(p(x, = sm sm sm — , (x, zjGii. (76) 

a a a 

Let 

A = I (logno)(^dxd2;, B = 2 [ ipoz^xixdz, loo = ( , ■ 
Jn Jn Jo Vx-^ + 1 

If A > and B > 0, then the solution of |7^[j-([73|) will blow up no later than T* = 
\18) 
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Proof We prove the theorem by contradiction. Suppose that the system (j74p has a 
globally smooth solution, li, V'z G C^{[0,oo); {{"^{Q)). Multiplying (p^ to the both sides of 
the ■i/'-equation and integrating over $7, we get 

- / AV't0zdxdz= / {u^)^(l),dyidz. (77) 
Jn Jn 

Note that = as long as the solution remains smooth. By integrating by parts and 
using the boundary condition on ^ and the property of (j) to eliminate the boundary terms, 
we have 

— / A^dxdz = / u^c/i^xdxdz. (78) 
Jn Jn 

Substituting </> to (j78p and using the boundary condition for ip, we obtain 

^z4'&x.dz = - n^^dxdz. (79) 



/ (logM)</>dxdz = 2 / Vz</>dxdz. (80) 
Jn Jn 



Jn 3 
Next, we multiply </> to ([8]) and integrate over J7. We obtain 

d. 

Combining ([79]) with ([80]) . we obtain 

d^ /■ 2 /■ 

/ (logn)(/)dxd2; = - / n^(^dxdz. (81) 

Integrating the above equation twice in time, we get 

/'(logn)(/)dxdz = -[ [ ( [ u'^(j)d^dz\dTds + A + Bt. (82) 
Jn 3 Jo Jq \Jn J 

Using (j82p . following the same argument as in Theorem 13 -H we can prove that the solution 
of the initial boundary value problem (j74p blows up no later than T*. □ 



4 Blow-up of the 3D model with partial viscosity 



In this section, we prove finite blow-up of the 3D model with partial viscosity. Specifically, 
we consider the following initial boundary value problem in a semi-infinite domain: 



-Alp 



UJ. 



(x, z) G = X (0, oo) 



The initial and boundary conditions are given as follows: 

^\dn\v = 0, {i^z + /3V') Ir = 0, 
^\dn\v = 0> {'^z + 1^) |r = 0, 
Lo\t=Q = ioo{yi,z), M|t=o = ^io(x, 2) > 0, 



(83) 



(84) 
(85) 
(86) 
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where F = {(x, z) € | x G ilx, -2 = 0}. 

Now we state the main result of this section. 

Theorem 4.1 Assume that uqIqq = 0, uo\q > 0, and UQ,ipQz £ H'^i^), £ H^(Q). 
Further, we assume that [3 G 6*00 as defined in Lemma \2.1\ and [3 > 7 = Let 

(/.(x, z) = e-"' sin — sin — , (x, z) e Q, (87) 
a a 

where a = satisfies < a < v^vr/a. Define 

A = / (log uojc/xlxdz, B = — ojQcfizd'Kdz, D 
Jq Jn 



dx _ fira^D'^B^ 



// A > 0, 5 > 0, and r* < (log 2) (v - aA) \ then the solution of model (8^ with 



im 



tial and boundary conditions 184^-^8^) will develop a finite time singularity before T* . 



Proof Again, we prove by contradiction. Assume that the 3D model with the initial 
boundary conditions ([8^ - ([86]) has a globally smooth solution, u, ipz S C^([0, 00); if^(J7)), 
u) G C^{[d,co)]H^{yt)). Multiplying to the both sides of the •i/'-equation and integrating 
over 0, we get 

- / A'i/'^zdxdz = / cj^^dxda:. (90) 
Jn Jn 

By integrating by parts and using boundary conditions ([8^ - ([85|) and the property of (p, we 
obtain 

/ 'i/'zA^dxdz — / '02</'2|2=odxdz — / -0^x01 2=odxd2; = / cjc/i^dxdz. (91) 
Jn Jn^ Jn^ Jn 



Substituting (p defined in (jSTj) into the above equation, we have 



0i(x)dx 

2=0 



a / V'2<Adxd2; = / ^(/i^dxdz - / aipz H 

. « J Jn Jn Jn^ V a y 

= j^iocPzdxdz + V'|2=o0i(x)dx, (92) 
where ^1 (x) = sin sin . Since /? > , we can choose 



a = ^ < (93) 



27r2 V2vr 

to eliminate the boundary term in (j92p . This gives rise to the following identity: 



■ipz(t><i:x.dz = 2 / i^02dxd2:. (94) 



n 2(f) -a2Jn 
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Next, we multiply 4)z to the both sides of the w-equation and integrate over Q, 

1 wt^zdxdz = / {u^^ _^(j)zAx.dz + V I Awc/i^dxdz. (95) 

Integrating by parts and using u\qq = 0, we obtain 

/ uj<j)zdxdz = - I u^(t)zzd-x.dz + v (- [ uJz(l)z\z=odx + [ ^^4>zz\z=odx + [ wA(/)^dxdz) 
Jn Jn \ Jn^ Jn^ Jn / 

= —o? I (f)d:>cdz — u ( — a^j / wc/i^dxdz + f {aojz + a^oj)\z=o(t)i{'x.)d-x 

Jn V « J Jn Jn^ 

= —c? j n^0dxdz — I — ^ oij [ w^^dxdz + / a{a — ■j)uj\z=o(f>i{'x.)d'x 

Jn V « J Jn Jn^ 

= —a^ / (pd-Kdz — u ( a^] [ wf/i^dxdz, (96) 

Jn V « / Jn 

where we have used a = 7 to eliminate the boundary term in the above estimates. Solving 

the above ordinary equation for (j^^dxdz gives 

/ uJ^zd:>^dz = e"^* f L^o'/'^dxdz - [ e'^^*"") ( [ u^^d^dz^ ds, (97) 
Jn Jn Jo \Jn J 

where A = Using the renormalized u-equation ([8|) , ([M|) and ([97|) , we obtain 

/ (logu)(/)dxdz = 2 / ^pz(pdyLdz (98) 
dt Jn Jn 

e"^* [ Loocpzdxdz + [ e"^^*"') ( [ u^^dxdz] ds 



2{l)'-a^\ Jn ' Jo 
Integrating the above equation in time, we get 

/ (logu)(f>dxdz = I (logt(o)</'dxdz — | (- — ^ ) (- / cJo'/'zdxdz) 

Jn Jn 2(^) - \ J \ Jn J 



2(f)' 

+ r7r^? I r e-^(^-") ( / nVdxdz ) drds. (99) 



2^2 rt rs 



2i^y -a^ Jo Jo 



Let To = then e"^* > ^ over the interval [0,ro]. Note that ^ (^-=x^) = e"^* > ^ 

A - 2 



for < t < To. This implies that M — > I for < t < Tq. Thus we have from (HH) that 



/"(log«)0dxdz > A + l-DBt + ^Da'^ [ [ ( [ u^cpdjcdz^ drds, (100) 
Jn 2 2 Jq Jq \Jq J 

for all t G [0, Tq]. Now we can follow exactly the same procedure as in the proof of Theorem 
3.1l to prove that the 3D model must develop a finite time blow-up before 



^•=(-12^; w 

Since T* < Tq, we conclude that the solution must blow up before T* . □ 
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Remark 4.1 We can also prove the finite time blow-up of the 3D model with partial vis- 
cosity in a bounded domain following a similar argument. We omit the analysis here. 



5 Blow-up of the 3D model with conservative boundary con- 
ditions 

In this section, we will consider boundary conditions for ip that will conserve energy. Under 
some additional condition, we can prove that the solution of the 3D model with conservative 
boundary conditions will also develop a finite time singularity. 

5.1 Blow-up in a semi-infinite domain 

Consider the following initial boundary value problem: 

I a"^* " ' (x,z)Gf^ = flxX(0,oo), (102) 

i^\dn\r = 0, i^zlr = 0, (103) 
V'|t=o = V'o(x, z), u\t=o = no(x, z) > 0, 

where fix = (0, a) x (0, a), and F = {(x, z) G | x € Ox, z = 0}. 

Theorem 5.1 Assume that uqIqq = , uqIq > 0, and uo,^oz G H^{0,). Let 

(/.(x,z) =e-""sin — sin— , (x,z)eO, (104) 
a a 

with a = ^, and 

A= (log no)</>dxdz, B = 2 Voz^dxdz, 
Jn Jn 

^(^) = . x2 / (V' - '0o)U=osm sm dx. 

2(f) -a2yn, a a 

If A > 0, B > and r{t) < ^ as long as u.,'4) remain regular, then the solution of II 10^) - 
hl 03\) will develop a finite time singularity in the norm. 

Proof We prove the theorem by contradiction. Assume that the initial boundary value 
problem has a globally smooth solution, u, ^2 G C-'^([0, 00); //^(l])). Multiplying to the 
both sides of the '(/'-equation and integrating over fi, we get 

- / AV't02dxdz= / (u2)^(/>2dxdz. (105) 

Note that u\z=o = since mo|z=o = 0. By integrating by parts and using the boundary 
condition of ip and the property of (j), we have 

- / Vzi'^'/'dxdz - / 'ipzt4>z\z=od-^dz - / Ax0U=odxdz = / u^^^^dxdz. (106) 
Jn Jn^ -/Hx Jn 
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Substituting (p defined in (jl04p into the above equation, we have 

(^2(^)'-a2^ ^^V'.0dxdz = a2^n2,/,dxdz-2(^)'^ V'tU=o</'i(x)dx, (107) 
where (pi (x) = sin sin . Finally we have 

^ [ V'.^dxdz = — / u^dxdz- '^[^^ [ Vd.=o</'i(x)dx. (108) 

Next, we multiply (j) to ^ and integrate over $7. We get 

[ (logtt)(/>dxdz = 2 [ ip^(t)d-x.dz. (109) 
Combining (jlOSp with (jl09p . we obtain 



dt2 



/■ 2q2 /■ 4('-)^ r 

/ (log^z)</.dxdz= — -2 / n^^dxdz-— / VtU=o0i(x)dx (110) 

Jo, 2 (^) - 2 (^) - 



Integrating the above equation in time and using the assumption that r(t) < ^, we get 



/ (log u)(j)Ax.Az 
Jn 



2a^ 



2 (f)^ - q2 Jo Jo \Jn 



t rs 



n^Mxdz ) drds + A + Bt 



2{K)'-a^Jo \Jn. 



{ip - -00)12=0 '/'I (x)dx I ds 



2[^'-a^Jo Jo \Jn 



> ^ / / f / nVdxdz ) drds + A + ]-Bt. (Ill 
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Using (imp , following the same argument as in the proof of Theorem I3.H we can prove 
that the solution of the initial boundary value problem of the 3D model blows up in a finite 
time. □ 



5.2 Blow-up in a bounded domain 

In this subsection, we will prove the finite time blow-up of the 3D model with a conservative 
boundary condition in a bounded domain. Specifically, we consider the following initial 
boundary value problem: 

j " ' (x,z)ef7 = QxX(0,6), (112) 

i^\dn\r = 0, ^z\r = 0, (113) 
V'|t=o = ^o(x, z), u\t=o = tio(x, z) > 0, 

where x = {xi,X2), fix = (0, a) x (0, a), F = {(x, z) E | x G ilx, z = or z = 6}. 
The main result is stated in the following theorem. 
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Theorem 5.2 Assume that uq\qq = , uo\q > 0, and uo,^poz G H'^{Q). Let 

(p{x,z) = sm sm , (x, zjGiZ, (114) 



with a = f , and 



A= {\oguo)4'd'Kdz, B = 2 ■i/'o^f/'dxdz, 
Jq Jq 

r{t) = ^^^2_^, y^j^- Vo)U=osm — sm— dx< -. 

7/ A > 0, B > Q and r{t) < ^ as long as Ujijj remain regular, then the solution of ill^} - 
illS]) will develop a finite time singularity in the norm. 

Proof We prove the theorem by contradiction. Assume that the initial boundary value 
problem has a globally smooth solution, «, V'z £ C^([0, oo); ff^(O)). Multiplying (f)^ to the 
both sides of the ^/'-equation and integrating over 57, we have 

- / AV't^^dxdz = / (u^)^(/)2dxdz. (115) 

Note that u\r = as long as the solution remains regular. By integrating by parts and 
using the boundary condition of ^jJ and the property of 4>, we get 

- / V^tA(/)dxdz+ / iptA^(p\lzld:>^dz = / ^^(/.^^dxdz. (116) 

Substituting (p to the above equation, we obtain 

I 2 f— ) — ) — / 'ipzCpdx.dz = [ u^(j)dx.dz 
\ \aJ J dt 

_P)'(e"^_e-"'') /■ VtU=o0i(x)dx, (117) 

where 6-1 (x) = sin sin . Thus we have 

1 /■ 2 /• f— )^ (e^^ e~'^^) I 

- / ^,(/,dxdz = — -2 / u^(t>d^dz- ^2 / V'tU=o0i(x)dx. (118) 

dUn 2{l) -a^Jn 2 (f ) - a2 Jf^x 

Next, we multiply 4> to ^ and integrate over We have 

/ (logn)(/.dxdz = 2 / V^^(/>dxdz. (119) 
dt Jn Jn 

Combining (jllSp with (jll9p . we obtain 



-2 / (logn)<^dxdz = -— 2 / nVdxdz- / ^i|,=o<Ai(x)dx. (120) 
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Integrating the above equation in time and using the assumption that r(t) < ^, we get 

/" (lognl^dxdz = f f ( f u'^M-Kdz] drds + A + Bt 

Jn 2(^?-a^ Jo Jo \Jn J 



7r'\2 f ^ab ^—ab\ ft ^ ^ 

{ip - -00) U=o'/'i(x)dx ) ds 



2(f) (e--e 



2(f) 7o \Ja 

> I I I I u^<l)dxdz] drds + A + Ist. (121 



2(f) -«^^o Jo \Jn 
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Using p21|) and following the same argument as in the proof of Theorem I3.H we can prove 
that the solution of the 3D model will develop a finite time singularity in the norm. □ 



5.3 Blow-up of the 3D model with other conservative boundary condi- 
tions 

The singularity analysis we present in the previous subsection can be generalized to study 
the finite time blow-up of the 3D model with the same boundary condition along the xi and 
X2 directions as in Section [5.21 but changing the Neumann boundary condition along the z- 
direction to a periodic boundary condition. The assumption on uq and -00 remains the same 
as in Section 15. 2[ In this case, we can prove the finite time blow-up of the corresponding 
initial boundary value problem with two minor modifications in the statement of the blow- 
up theorem. The first change is to replace (j) by the following definition: 

(l)(x,z) = - — — ^ ^ sin — sin— , (x, z) G 17 = (0, a) x (0, a) x (0, 6), (122) 

2 a a 

with a = ^. The second change is to modify the definition of r(t) as follows: 

2a(l — e~"*) /" , , vrxi Trx2 B 

2 (f) -a^ Jn^ a a 2 

where A and B are the same as in Theorem [521 If A > 0, 5 > and r{t) < f as long 
as n, ip remain regular, then we can prove that the solution of the corresponding initial 
boundary value problem will develop a finite time singularity in the norm. 

The same singularity analysis can be applied to study the finite time blow-up of the 
3D model with the same boundary condition along the xi and X2 directions as in Section 
15. 2[ but changing the Neumann boundary condition along the z-direction to the Dirichlet 
boundary condition. The assumption on uq and ipQ remains the same as in Section 15.21 In 
this case, we can prove the finite time blow-up of the corresponding initial boundary value 
problem with two minor modifications in the statement of the blow-up theorem. The first 
change is to replace (p by the following definition: 

^a{z-b) _|_ ^-a{z-b) ^ ^ 

(t>U,z) = sin — ^sin — (x, z) G = (0, a) x (0, a) x (0, 6), (123) 

2 a a 
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with a = ^. The second change is to modify the definition of r(t) as follows: 

a(e"^-e-"*) f, , , irxi 7rx2 , B 

^(*) = rT2 TT / - i^Oz)\z=o sm sm dx < — , 

2 (f ) -a^ Jn^ a a 2 

where j4 and B are the same as in Theorem O If A > 0, 5 > and r(t) < f as long 
as u, ip remain regular, then we can prove that the solution of the corresponding initial 
boundary value problem will develop a finite time singularity in the norm. 

Remark 5.1 All the results in this section can be generalized to a cylindrical domain 0, in 
high dimension space M^, with Q, = |(x, z)! x e Ox C M^~^,2: E [a, 6] C M}. In this case, 
the weight function 0(x, z) is chosen to be the product of two functions: 

(/)(x,z) =0i(x)r/(z). (124) 

Here the eig en- function, r]{z), is the same in the Eulerian coordinate in the previous sections. 
The eigen-function, </'i(x), defined in the x space, is chosen to be the first eig en- function of 
the following eigenvalue problem: 

-Ax0i = A(^i, (125) 
^ilan^ = 0, (126) 

with A > 0, where Ax is the N — 1 dimensional Laplace operator. Ax = + • • • + 



dx^ dx^ 

6 Global regularity of some initial boundary value problems 

In this section, we will prove the global regularity of the 3D model for a special class of 
initial and boundary conditions. To simplify the presentation of our analysis, we use and 
ipz as our new variables. We will define v = tf^z and still use u to stand for u^. Then the 3D 
model now has the form: 

I " ^™ , {^,z)en = {0,S)x{0,6)x{0,6). (127) 

-Avt = Uzz 

We choose the following boundary condition for v. 

v\dn = -4, (128) 

and denote v\t=o = vo{x, z) and u\t=o = uoix, z) > 0. 

In our regularity analysis, we need to use the following Sobolev inequality |12j : 

Lemma 6.1 For all s G Z+, there exists Cg > 0, such that, for all u,v e L°° Ci {^^) , 

\ 1/2 

J2 ||5"(nt>) -a^tt-ull^a I <Cs{\\u\\lo-\\v\\h- + \\Vv\\l--\\u\\hs-i). (129) 

,0<|a|<s / 
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Now we state the main result of this section. 

Theorem 6.1 Assume that uo,vo G H^{0,) with s > 4, uq\q^ = 0, and vq < —4 over 
then the norm of u, v remains bounded for all time as long as the following holds 

6{4Cs + 1) {\\vo\\h^ + CsWuoWh^) < 1, (130) 

where Cg is an interpolation constant. Moreover, we have \\u\\loc < ||tio||L°°e~'^*, < 
||uo||_H's(f7)e~'^* and \\v\\H^{n) ^ C for some constant C which depends on uq,vo and s only. 

Proof First, we note that (jl30p implies that 2Cs(5||uo||h'' < ^- By using an argument 
similar to the local well-posedness analysis, we can show that there exists Tq > such that 
llfxlliys and ll'^ll/f" are bounded, v < —2 on Q, and 2Cs(5||^;||/fs < 1 for < t < Tq. 

Let [0, T) be the largest time interval on which and are bounded, and both 

of the following inequalities hold: 

V < -2 over n, 2Cs6\\v\\h'^ < 1. 

We will show that T = oo. 

For a = (01,02,03) with oj > {j = 1, 2, 3) and |a| < s, we have for < t < T that 

^{d''u,d'^u) = 8{d''{uv),d'^u) 

= 8 {d°'u ■ V, a"n) + 8 (d^iuv) - d^'u ■ v, d^'u) 
a^npvdxdz + 8 {d'^iuv) - d'^u ■ v, 



< -16 / |a"n|2dxdz + 8||a"(tiw) z;||j^2||9"n||i2. (131) 



Using Lemma |6. 11 we get 
d 



iiii||_f/» < — 8||u||hs + (||ti||L°o + IIVvIIloo (132) 
at 



Since u\qq^ = "UolsfZ = 0, we obtain 



u{'x.,z,t) = / dziu{ii., z' ,t)dz' 

rz pxi i'X2 



d^i d^t dziu{x'i,x'2, z' , t)dxi dx2 dz' 
10 JO JO ' ' 

< 6'^^^\\dc,^d,,d,u\\L'^<S\\u\\Hs, (133) 



since s > 4. Notice that Vxi\z=o = 0; so we have 



v^^:,'dz' < I \vx,z'\dz' < 6\\v.j;^:,\\l^. (134) 
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Similarly, since Vz\xi=o = 0, we have 

Vz = I v^: ,<^xi< / \v^i ,\'^xi < 5\\vr,^z\\L^ ■ (135) 
Jo Jo 

Combining (jl34p with (|135p . we get 

\\Vv\\l^ < 6max{\\vxiz\\L°°)- (136) 

2=1,2 

Since s>4> 2 + 3/2 by our assumption, we obtain by using the Sobolev embedding 
theorem [12] that 



ll^'xizlli'^ < Cs\\Vx^z\\hs-2 < Cs\\v\\h'^- (137) 
It follows from and (fW]) that 

IIVvIIloc < C^^ll-yll//.. (138) 
Combine (fT32]) - (fT33]) with (fT38D . we obtain 

< (-8 + 2C,,5||?;||hO I^IIh- (139) 
Since 2Cs5||?;||//s < 1 for t < T by the assumption of T, we have for t <T that 

\\u\\h^ < \\uo\\H^e~'^^. (140) 

Note that 

^{d''v,d°v) = 2{d'^vt,d°v) <2\\d''vt\\L2\\d'^v\\L2. (141) 
Recall that Avt = Uzz- We can easily generalize the proof of Lemma |2. II to show that 

ll^^t||_H"»(f7 < Cs\\uzz\\H=-^{n) < (^8 11^11/^3(0). (142) 

Using p42p . we get 

^ll^^ll^s < 2||'t;j||//s||i;||^3 < 2C5||n||//3||'t;||//s. (143) 
Substituting (jl40p to the above equations, we get for t < T that 

< CsWuWh" < C^||no||/f»e"'^*. (144) 

at 

Integrating the above inequality in time, we obtain the estimate of over [0, T): 

WvWh" < llwollfl'^ + C'sll^'ollfl'^ / e~'^'^ds 

Jo 

C 

< \\vo\\h- + ^\\uo\\h^ < WvoWh" + CsWuqWh"- (145) 
24 



Since v\Qn = —4, we can use the the same argument as in the proof of (jl33p to show that 

\v + M< 5\\v\\hs < 5 {\\vo\\hs + CsWuqWhs) , (146) 

where we have used (jl45p . Now we have for t < T that 

V < -4^ + 5 {WvqWh" + C^lhoIkO , 
2Cs5||v||h= < 2Cs5 (||folk- + Cs\\uo\\h'^) ■ 

By our assumption on the initial data, we have 

{ACs + 1) (WvoWhs + CsWuoWhs) 5 < 1. (147) 

Therefore, we have proved that if 

v<-2onn and 2Cs6\\v\\h^ < 1, 0<t<T, (148) 

then we actuaUy have 

V < -3 on 17 and 2CsS\\v\\h- < ^, 0<t <T. (149) 

This implies that we can extend the time interval beyond T so that (jl48p is still valid. This 
contradicts the assumption that [0, T) is the largest time interval on which ()148p is valid. 
This contradiction shows that T can not be a finite number, i.e. (I148P is true for all time. 
This implies that ||ii||_f/s(o) and ^-re bounded for all time. Moreover, we have shown 

that \\u\\l°° < ||'Uo||L°°e"'^*, < \\uo\\H={n)e~'^^ and < ||?^olk« + CsWuoWh" ■ 

□ 



Appendix A. Proof of Lemma 12.11 

Proof of Lemma 12.11 We present the proof for the case of a = vr. The case of a 7^ vr 
can be proved similarly. First, we perform the sine transform along xi and X2 directions to 
both sides of ([20]). We have 

\k\^v{k,z)-v,,{k,z) = f{k,z), (A-1) 



where k = (fci, ^2); 1^1 = \/ki + ^1 , and the sine transform of v is defined as follows: 



2X2 



vr 



v{k,z)=[ — ] / / v{xi,X2,z)sm{kiXi)sm{k2X2)dxidx2- (A-2) 



/o JO 

Applying the sine transform to the boundary condition gives 



{v,{k,z)+f3v{k,z))\,=o = 0. (A-3) 
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The second order ODE (|A^ can be solved analytically. The general solution is given by 



v{k, z) 



\k\ 



fe-\''\''dz' + Ci{k)] + 



\k\z / Y 



fe^'^^^'dz +C2{k)] . (A-4) 



1^1 V2 

The boundary condition ()A-3|) and the constraint that v G L^(Q) determine the constants 
Ci and C2 uniquely as follows: 

f{k,z)e-\'\''dz', C2(A;) = 



Ci{k) = ■ 

Let xi^) be the characteristic function 



Ci{k). 



0, X < 0, 

1, a;>0. 



(A-5) 



(A-6) 



Then v{k, z) has the following integral representation (note that P ^ \k\hy our assumption): 



1 /'°° 

v{k,z) = -^J^ f{Kz)e-\^\^^'-^^^(^' 



\k\+l5 



\mk\-f3) Jo 



)^^,' _ ,)dz' + ^ m z)e-\'^\^^-^'h{z - z')dz' 
f{k,z)e-\^\^'+^''^dz' 



1 



-— f{k,z)Ki{z'-z)dz' + — I f{k,z)Ki{z-z')dz' 



2\k\ 

J{k,z)K2{z + z')dz' 



(A-7) 



\mk\-p) Jo 

where Ki{z) = e~l'^l^x('2), K2{z) = e~^^^^. Using Young's inequality (see e.g. page 232 of 
we obtain: 

\k\+P 



\vik,-)\\mo,oo) < ^ (2||iri||ii[o,oo) + 2 



1 1 ^2 1 1 Li [0,00) ll/(fcr)llL2[0,oo) 



< 



where M = max ( 1 + 

fel,fc2>0 



1 

1^ 
\k\+P 



1 + 



\k\ + f3 



M 



II/(^,-)IIl2[o,oo)<^II/(A;,-)IIl2[o,oo), (A-8) 



p 



J 



< 00 since P ^ \k\ for any A; G Z'^ by our assumption. 



Next, we estimate Vz{k,z). Differentiating ()A-4p with respect to z, we get 

v,{k,z) = — / f{k,z)e-\^\^^'-'^dz' - - f{k,z)e-\^\^'-''Uz 



2 _ 
\k\+P 



f{k,z)e-\^\^^+''Uz'. 



(A-9) 



\k\-P. 

Following the same procedure as in our estimate for v{k,z), we obtain a similar estimate 
for Vz{k, z): 



I^^(^>-)IIl2[0,oo) 



< 



1 

W\ 



1 + 



\k\ + p 



p 



M. 



II/(^>-)|Il2[o,oo) < j^ll/(^r)|lL2[0,oo) 



(A-10) 
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Let a = (01,02,03) with aj > (j = 1,2,3). We will prove ||9°?;||22 < W^Wff^^^^^^ for 
all \a\ > 2. We will prove this using an induction argument on 03. First, we establish this 
estimate for 03 = and 01+02 ^ 2. Below we use the case of oi > 1 and 02 > 1 as an 
example to illustrate the main idea. By using the Parseval equality and ()A-8p . we obtain 

00 00 „oo 
tl = lfc2=l •'^ 



k 



00 00 



= E E^i"^Mi^(^'-)iii^[o,oo) 

ki=l k2=l 

00 00 

< E E^'^i"^^M^i~'ii/(^'-)iii^[o,oo) 

fci = l k2=l 

00 00 

< ^||C"'^"^-7(fc,-)|li2[o,oo) 

fci=l k2=l 

= M2||a:-ia,"-Vlli2(n) <m2||/||^i„i_,(^^. (A-ii) 

Similarly, we can prove (|A-lip for 03 = and oi + 02 > 2 by distributing the appropriate 
order of derivatives to xi and/or X2 direction. 

Using (1A-10|1 and following the same procedure as in the proof of (lA-lip , we can prove 
()A-lip for the case of 03 = 1 and oi + 02 > 1. Finally, using (|A-ip and differentiating (1A-1|) 
with respect to z as many times as needed, we can prove 

\\d''v\\h^n)<Ca\\f\\W.^^y (A-12) 

for all 03 > 2 and oi + 02 > by using an induction argument and (lA-lip for 03 = and 
03 = 1. Using (1A-12P and (lA-7p . we obtain 

\MH-{n) < Cs\\f\\H--2{n), (A-13) 

for all s > 2, where Cs is a constant depending only on s. The uniqueness of the solution 
follows from the solution formula ()A-4p and ()A-5P . This completes the proof of Lemma 12.11 

□ 
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